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I. INTRODUCTION. 

Canonical Perturbation Theory (CPT) was historically the first perturbation theory in 
quantum mechanics 1 . It iteratively constructs a unitary transformation of the Hamiltonian 
matrix into a block-diagonal form. The resulting effective Hamiltonian 2,3 then allows for a 
solution of the eigenvalue problem. 

Nowadays, Van Vleck 4 CPT in Primas’ superoperatorial formalism 5 is widely used in 
quantum chemistry and molecular spectroscopy 6,7 . Contemporary precise models demand 
time-consuming computations up to very high perturbative orders 8 . This is why new efficient 
algorithms are important. 

We present here a noniterative construction of the block-diagonalizing transformation 
for CPT. Our approach follows the recent application of Kato perturbation expansion to 
classical mechanics 9 . It uses the remarkable analogies between mathematical formalisms of 
perturbation expansions in classical and quantum mechanics. 

The method is based on a Kato series 10 for the the Laurent coefficients of the resolvent 
of the Liouville superoperator and the Dyson’s ordered exponentials. We demonstrate a 
regular structure in perturbation series and derive an explicit expression for the generator 
of a block-diagonalizing transformation in any perturbation order: 

G = iS A H 1 . 


Here, Sh is the partial pseudo-inverse of the perturbed Liouville superoperator. Unitary 
intertwining of perturbed and unperturbed averaging superprojectors allows for a description 
of ambiguities in the generator and block-diagonalized Hamiltonian. The corresponding 
perturbative algorithms are sufficiently efficient even for high-order computations. 

Note that all calculations here are only formal, in a sense that neither a discussion of 
power series convergence nor conditions for the existence of constructed superoperators are 
present. 

The downloadable “Supplementary data hies” 11 contain demonstrations and large formu¬ 
lae. The demonstrations use the freeware computer algebra system FORM 12 . 
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II. BASIC PERTURBATION SUPEROPERATORS. 

Consider the Hamiltonian H that differs from an exactly solvable one by a perturbation: 

H = H 0 + ctH x . 

We will discuss here only time-independent Hamiltonians that have a discrete spectrum. 

The goal of canonical perturbation theory is to transform the perturbed Hamiltonian 
into a simpler block-diagonal operator by a near-identity unitary transformation. Such 
transformations may be handled conveniently using superoperator formalism 5,13 . Here we 
will outline only some basics of it. 

For any G, one may introduce a quantum Liouvillc superoperator acting on the operator 
space: 

A = [G, .], C e F=[G,F]- 

We will use the bold calligraphic-style letters here to denote the superoperators. Due to the 
Jacobi identity, the Liouvillian is the “derivation” of a commutator: 

C 6 [F,H} = \CaF,H] + [F,C 6 H]. 

The C of the discrete spectrum system admits the decomposition of its domain into the 
direct sum 91© 91, where 91 or 91 are the kernel and the range space, respectively. In order to 
realize this decomposition, Primas 5 introduced the following averaging (block-diagonalizing) 
superoperator: 

/*+oo , a /*+oo 

V fr F = lim A / e~ xt e^° Fe~^°dt = lira A / dte~ xt e* tCft oF. (1) 

H ° a^+o J 0 a^+o J 0 v ; 

-2 

This superoperator is the projector 'Pfio - ^Ho- It projects the operator F onto the kernel 
space 91 of block-diagonal operators commuting with H 0 . 

The complementary projector 1 — "Pgj Q extracts the non-commuting with the H 0 (off- 
diagonal) part of the F. It projects on the range space 91 where the inverse of Cjg Q exists. 
This inverse is the integrating superoperator <Sh 0 , which is also called the solution of the 
homological equation, tilde operation, zero-mean antiderivative, Friedrichs T operation, £ 
operator, division operation, etc. Its invariant definition by Primas 5 is as follows: 

p OO . 

A„ = -WW o *e-"e*‘ £ *°( 1-A 0 )- P) 
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It is easy to check that: 


^Ho^Ho ~ ^Ho^Hq ~ 1 ^H 0 i 

= °- 

The superoperators {C-p x , F g 0 - S h 0 ) are the building blocks of the canonical perturbation 
series 5 . Quantum mechanics uses the following methods for their computation: 

• The Energy Representation of an unperturbed Hamiltonian provides the eigenbasis of 
the Liouville superoperator 13 . The action of the basic superoperators on any operator 
F can be directly computed in this representation as follows: 

F 

J- m.r 


'Ph 0 F= ^ F mn\m) 


n\ 




E 


Ffn E n 


| m) (n |. 


Em=E n E m ^E. 

These well-known expressions are used frequently as the definitions of Fp Q and S. 

The method of Contact Transformations 14 (quantum Birkhoff-Gustavson normal 
form 15,16 ) applies to the perturbed harmonic oscillator 

d 

k ( ~2 


H = (pi + ® + a#i(p, q). 


k =1 


It uses the representation of ladder operators: 


Qk = \l^{a'k + a k), 


Pk 


= (4 - a k ), [a k , aj] = S M . 


The action of basic superoperators on a Wick-ordered polynomial operator F(p, q) = 
F(m, n) a) m a n is as follows: 

T>i h F= E S 6 F= V - F( “’ S) 

(tjj,rh—ri)=0 n):/0 


h(uj, rh — n) 


a^a H . 


III. THE SUPERRESOLVENT 


Consider the resolvent of the Liouville superoperator: 

1 


&h(z) = 


-' C’H — Z 


( 3 ) 


This superoperator-valued function of the complex variable z is the Laplace transform of 
the evolution superoperator 
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r+oo 

& H (z) = - dte~ zt eP CH . 

Jo 

Resolvent singularities are the eigenvalues of the Liouvillian. Typically, its spectrum 
a(Cfi) = ff(H) — cr(H ) has a richer structure than that of the corresponding Hamiltonian 17 . 

Let us begin with the simple case of the pure point spectrum of a quantum system 
with isolated energy levels. The singularities of corresponding Liouvillian superresolvent are 
located at E k — E m and are separable from the origin. 

The existence of E^ and <Smeans that the unperturbed superresolvent has a simple 
pole in 0. The averaging superoperator is its residue in this pole: 

'Ph 0 = ~ R es z=o &h 0 i 

while the integrating superoperator S is its holomorphic part: 

S h 0 = 1 1 }™^h 0 ( z )( 1 -'Ph 0 )- 

Therefore, the Liouvillian sup erresol vent combines both basic perturbation superoperators 9 . 
This allows us to apply the powerful formalism of complex analysis to perturbation theory. 
It is well known 10 that, due to the Hilbert resolvent identity 

&h{z l) -&h( z 2) = (-1 - Z2)&a(zi)£ii(z2), 


the Laurent series for the unperturbed resolvent has the form of a geometric progression 

ra+l +00 


^H 0 ( Z ) ~ H 0 


E 

n =0 


h -S 


Ho 


E(t 


n=0 


M re—1 


J^H 0 Z 


Here, we have denoted 7 'Z^ 0 = ~'Ph 0 and 'R^l = $h 0 - 

The perturbed resolvent may be more singular. The Laurent series for a general resolvent 
with an isolated singularity at the origin has the following form 10 : 


&h ( z ) 






Here, is the eigen-nilpotent superoperator, which does not have an unperturbed ana¬ 
logue. 
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A. Kato series 


If the perturbation H\ is relatively bounded with respect to H 0 , then the Rayleigh- 
Schrodinger series converges 10 and the perturbed superresolvent can be expanded into the 
Neumann series as follows: 

^H 0 +aHS Z ) = ^H 0 ~ (^j ^hj~' hJ-' H 0 + ' ' ' 

= 5Z ( _0 +) ^h 0 ^ z ){^hJ^h 0 { z )) ■ 

n =0 ' ' 

The integration around a small contour results in the Kato series 10 for the “perturbed 
averaging superoperator”: 


i>H = 


2ni 


l*l=* 


^H{z) dz 


2ni 




n =0 
oo 


=r> J \ z \= e 


a h ’ A ° [z] 


Cpf^Mf_j Q {z) ) dz 


—E 

2ni ^ 


n=0 


l 

- a n 


\z\=e 


£ (?) “as 1 *- 1 

\m=0 ^ ' 


k =0 


j- \ k 

h\ 


4E < 


the “perturbed integrating superoperator”, and the “perturbed eigen-nilpotent’ 


1 H 


z ^H(z)d z i 'D ft — 2 n 


h 


z @H{z) dz - 


2l\h J\z\=e ' ' ^ J\z\ = 

Only coefficients of z _1 in these expansions will contribute to the result: 

( \ 

rv~i 

71+1 _ n \ ^ /+■» (P n +1) 




a 


,(Pn) 


,(P2) 


(Pi) 


E *!T'**** •■■****** 



n =0 

E Pj= n 
\ pd° 



n times 

y 1 




( 






\ 


= E(-1)" a” 

E 




■ ■ ■ ,X -H 0 t-H 1 



n =0 

EPi=n+l 

V pj>° 



n 

times 


/ 



( 






\ 


= E(-i rv 

E 




1 -^ 

JS<£t? 

• 1 



n= 1 

EPr”- 1 

V pj>° 



n 

times 


/ 


Because the superoperator consists of r p 3 superoperators Sp o , the summation in 

the above expressions should be done by all possible placements of n (or n + 1, or n — 1) 


dz , 


(4) 
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<S f l() in n + 1 sets. There are CV ln terms of order a n in Vfy and C 2 n+1 terms in For the 
Hermitian G, the operators S lf G and are anti-Hermitian. 

The properties of unperturbed superoperators can be extended to their analytic contin¬ 
uations as follows: 

VfjH = H , S = 1 — 'Pfj., Cfi'Pfi = VfjCfi = T>ff, etc. 

For the details, see our work 9 on an application of the Kato series to classical mechanics. 

To avoid misunderstanding, it should be noted that 'P^F will not be commuting with 
the perturbed Hamiltonian. This is because C'h'Ph = T>ff 0 in general. Actually, the 
“perturbed superprojector” Vfy projects onto some analytic continuation of the algebra 
of integrals of the unperturbed Hamiltonian. This may not coincide, in general, with the 
algebra of integrals of the perturbed system, because of a destruction of symmetries. In other 
words, the zero eigenvalue of the Liouville superoperator may be split by the perturbation. 

B. Canonical properties of the Liouvillian superresolvent 

Because C is a derivative, there exists an integration by parts formula for its superresolvent 9 . 
For any z i, Z 2 , and Z 3 outside of the spectrum of and any operator F, the following 
holds true: 

h( z i) a (z 2 )f— & (z 2 )f^ h( z z) + &h( z i)£'f^'h( z 3) 

= ( z\ — Z 2 — Z 3 ) & fi(zi) & fj(z 3 ). (5) 

This results from the application of the identical superoperator £% h( z i)(^h — Z\) = 1 to 
the commutator ^(^ 3 ) and the expansion of the Jacobi identity. 

Consider the derivative of the Liouvillian superresolvent with respect to the perturbation 

^a k [z' i = ~at k [z)t^a k (z). 

Substituting z \ = z 3 = z , and F — Hi into the canonical superresolvent identity (5) yields 

3&h( z ) = t-'fc ^(22)^1^ h( Z ) ~~ ^h( Z ) ^(z 2 )Ai — Z 2^h( Z )£'£ 6 (z 2 )H 1 ^h( Z )- (®) 

Look at the coefficient of * 9 i n the Laurent series of the above expression: 

~Q^^h( Z ) — ^’h( Z ) — H^ z )^vfjHi ^h( Z )’ 
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Proceeding similarly for coefficients of z 2 n (n > 1) in (6), we obtain 
^'h( z )£"P 6 h 1 = ^h( z ) ~ 

= h( Z ) ~ ^^H( Z )^-'±> n ^+ 1 H 1 ^H( Z )- 


This allows us to rewrite the expression for the superresolvent derivative as 


d ( y^v( Z ) —^h( Z )^S 6 Hi h( z ) 

~ fjHi^ h ( Z ) 2 + (ft) fc-DjjHi&H 


[zf- 


h I h( z ) A + • 


Actually, this is a power series because T>p, = O(o n ). 

It follows from the Hilbert resolvent identity that 

r) n -hi 

= n\ a 2 (z). 


dz 


Finally, we obtain 


^h( z ) =®h£s 6 h 1 

+ l(i\ C 


2 \ h i ' D ^ Hl 


hJ^h( z ) 

- -C - 

h 

dfflfi{z) 

Hl dz 

SW H (z) 

1 (*' 

A S^ 6 (z) 

d z 2 

6 U 

J L "b\H i Q ~3 


+ ... 


(7) 


The derivative of superprojector {d/d a)V ^ can be obtained as the residue of this expression 
at z = 0. In our case of an isolated point spectrum, the superresolvent is a meromorphic 
function, and therefore the residue of any of its derivatives with respect to z vanishes iden¬ 
tically. As a result, the superprojector 'P^ transforms canonically under perturbation: 


C s k H^H- 


This identity holds to all perturbation orders. 


( 8 ) 


IV. AN EXPLICIT EXPRESSION FOR A GENERATOR 
A. Ordered exponentials 

Ordered exponentials were introduced in quantum field theory by Dyson 18 in 1949. They 
parametrize an a-dependent unitary transformation U{a) 

\il>(0)) = U${a)), H = U~ l H U , 
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and its inverse U' = U 1 using direct and inverse Hamiltonian flows in “time” a with some 
Hermitian generator G(ot): 

r\ r\ 

-?-U = iG(a)U, ^-U ] = -iU ] G. 
da da 

In classical mechanics the ordered exponentials are known as Lie-Deprit transforms 19 . Their 
applications to quantum canonical perturbation theory are discussed by Ali 20 and Scherer 21 . 
We adopt the notation by Suzuki 22 for the direct and inverse ordered exponentials: 


U = exp + ( i / Gds) : W = exp_(- 


Gds). 


The integral signs here have only notational meaning as a reminiscence of the Dyson series. 
We will mostly use the superoperatorial ordered exponentials that satisfy the equations: 

Au(a) = a C 6(a) 0, A V(a) = i VC 6 . (9) 

Analogous to the Hausdorff identity, they are factorizable as follows 23 : 

POL POL POL 

U(a)F =exp + ( i / C^ds) F = exp + ( i / G ds) F exp _(—i / G ds ), 

Jo Jo Jo 

POL POL POL 

M(a)F =exp_(—* / Cq ds)F = exp _(—i / G ds) F exp + ( i / Gds). 

Jo Jo Jo 

Therefore, the transformed Hamiltonian can be written as H = exp _(—i / 0 ° C^ds) H. 

In quantum physics the ordered exponentials are usually computed by the Dyson series 18 . 
But canonical perturbation theory uses their representation as a power series: 

POL OO POL OO 

exp + ( i / Cg ds) — Qi n U n , exp _(—i / Cq ds) — a n M n ■ 


n =0 


n =0 


Substituting these series and G(a) = i n to the equations (9), Deprit obtained 

the relations for the coefficients: 


U r 


. n —1 
l 

n 


tX c 


Gn. — k 


Ui 


Mr 


k =0 


n —1 


- E V|X 

n < ^ 


G n . — k — 1 ^ 


( 10 ) 


k =0 
19 


and developed the “triangular” computational algorithm 19 . 
Iterations of (10) yield non-recursive formulas 24 : 


U, 


V, 


E 


C;, 

2. Axn — m-y — 1 — rn 2 — ^ 


n>mi>m 2 >--->m r 


n 


E h> 


c, 

•\r+l —1 


mi 


'Gn 


G rni — 1 
m r 


2 ~ 1 Gn — m^—1 


n>mi>m2>->mr 


m r 




n 
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Here, the sum runs over all sets of integers (mi,..., m r ), satisfying n > m 1 > ■ ■ ■ > m r > 0. 
In the first orders: 


Q 2 , -2 


ex P + ( i J o C G ds ) = l + iaC G 0 ~ y( £ G 0 ~ iC G^ 


a 3 , 2.3 


? g (A? 0 * t-'Go^-'Gi 2f Cfj Q 2£q 2 ) + 0(a ), 


I * A A 2 

exp _(—i C 6 ds ) =1 - z a£ 6o - y(^G 0 + *£&) 


Ct A 3 /V /V /V /V /S 

* y (-^Go + 2z CqJLq x + z CqJLq^ — 2i^) + 0(a 3 ). 


B. Intertwining transformation 


The canonical identity (8) means that the ordered exponential with the generator z <S k H \: 

pot 

IV’(O)) = exp + ( z / iSfjHids) \ip(a)), 


H = exp_(—z 


■'z «S 


ds) H, 


unitarily connects the perturbed and unperturbed superprojectors: 

/ ct /»q; 

AyydsjP^exp-H JI £is 

Since this transformation intertwines the superprojectors 


ds). 


exp_(—z / ^ ds) = V ko exp_( 


'iS A H i 


ds), 


( 11 ) 


and P if II = id, the transformed Hamiltonian becomes block-diagonal: 


H = exp_(—z 


'iSf.Hx 


ds) V k H = V ko exp_(—z / 


ds) H = V k H. 


JO JO 

We have thus shown that the ordered exponential with the generator 




G = iS fl H 1 =i £(-l) 


n+1 


a 


n=0 


\ 


e **%,**- 

EPi=n+l V v '"" 




n times 


( 12 ) 


7 


\ *y>° 

formally block-diagonalizes the Hamiltonian in all orders in a. 

The above formula can be extended to systems with degenerate energy levels, or to 
unbounded perturbation in the form of an asymptotic series. It may be shown 9 that an ex¬ 
pression for G truncated at order N results in a block-diagonal Hamiltonian up to 0(a w+1 ). 
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For clarity, we here use formal “analytic” expressions. But for degenerate systems, it must 
always be remembered that these expressions must be converted straightforwardly into trun¬ 
cated sums, and all equalities hold up to 0 (qA +1 ). 

In the first perturbative orders, the generator 



S C fti VC H S VC 6 i SC Ai S VC k S C ki S-\-S C ki VC ki V 

+vt ll s ! 'c lli v + re. re.. s') it + O(o 3 ) 


yields the following block-diagonal effective Hamiltonian up to « 4 : 



Here we omitted the indices Ho for unperturbed operators *P ko and S ko - These expressions 
are more regular than corresponding formulas 6,7,25 for standard CPT methods. 


C. General form of the generator 

Knowing that V ko and V k are unitary connected and that this transformation block- 
diagonalizes the Hamiltonian, the procedure of canonical perturbation theory can be viewed 
as the construction of unitary transformation, which connects unperturbed and perturbed 
superprojectors. 

Let us now determine the general form of the generator of such a transformation. The 
perturbed V k should satisfy the differential equation 



Application of this expression to the Hamiltonian H and the identities 
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yield 

(l-V il )C il G = Hl-V lI )H 1 . 

To solve this equation, it is sufficient to apply the Sp superoperator. Therefore, the general 
form of the generator of intertwining unitary transformation is 


G = iSpH 1 + F A F, 


(14) 


where F(a) may be any Hermitian operator. This is the main result of this work. This 
formula provides a non-recursive expression for the generator of the block-diagonalizing 
transformation and defines its ambiguity. 

The choice of FpF is the uniqueness condition 26 . It is natural to choose F = 0 or 
FpG = 0. This is not equal to the “off-diagonal” condition Fp Q Gp = 0 traditionally used 
in canonical perturbation theory 5 . Because C^^pFp = FpCp, _ p, the superprojector F u 
is itself insensitive to this choice. 

We can conclude that the generators of block-diagonalizing transformations can differ by 
a function belonging to a continuation of the algebra of integrals of unperturbed system. 

Effective Hamiltonians that were block-diagonalized using different uniqueness conditions 
F\ and F 2 are connected by unitary transformation: 


U21 


pot. 

exp_(-d / Cq 2 ds) 
Jo 


pot 

exp + (i Cp i ds). 

Jo 


Let us find its generator: 


^-U 2 i = i exp _(-i 


**) (^Gi ) eX P+( * / ^Gi^ S ) 


— i exp_(—i / C d2 ds)Cp,,p i _p 2) exp + (i C 6l ds) 


= i C 


exp_(-»/“ ds)V A (A-F a ) 


u 


21 - 


Because of the intertwining relation (11), the generator of U 2 i is always block-diagonal: 

pot pot 

G21 — exp_(— i / £q 2 ds) F g(Fi — Ff) = F i j Q exp_(—i / C, is p 2 ds) (F\ — Ff). 

Jo Jo 

This is the well-known result that the effective Hamiltonian is determined up to block- 
diagonal transformation 2 . Such ambiguity does not affect the perturbed eigenvalues. 
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V. COMPUTATIONAL ASPECTS 

A major difference between this and standard canonical perturbation algorithms by Van 
Vleck 4 , Primas 5 , etc., is the explicit non-recursive formulas. Traditionally, perturbation 
computations block-diagonalize the Hamiltonian order by order. In contrast, we compute 
the generator G = iS 1 up to the desired order directly. Then, the ordered exponential 
block-diagonalizes the Hamiltonian. 


A. An explicit algorithm for the generator 

Let us introduce the superoperators 


£ m 

Z„ = 


(- 1 ) 


n+1 


£ 


n 


(Pn + 1) 

H 0 


(Pn) 


P 1 + ... +p n+ i=m 
Pj >0 


^hJ^Ho 


j^ V2) t - 'TZ' Pl ' > 

H 0 L 'H 1 H 0 ■ 


n times 

Here, the sum runs over all possible placements of m superoperators Sh 0 in n + 1 sets. The 

m 

Ho 


/s Q /\ /s 777 /■> 777 

simplest such expressions are Z {) = 'P Ho , and Z 0 = —S H for m > 0. We have already 


'0 r Hq 

discussed their computation in section II 
The following recursive relation 


~ m v—-v a m—k ~ K, 


k =0 


leads to an efficient algorithm for the generator G = iSnH\. We consider the square table: 


F 0 ° 

^o 1 

*o ■ 

t?N 

• • ^0 

pW+i 

HO 


Fl 

F? • 

.. Ff 

pN+1 


Fh -1 

F 2 n . i • 

riv 

• • r JV- 1 

P N +1 

r N- 1 


zpN+l 

Here the first row is 

F 0 ° = Z% = Vh 0 Hi, F 0 m (x) = ZqHi = ~S m Ho H i, 
and each subsequent row is generated from the previous row according to the rule 

m 

fpm _ \ ' fpk 

^n+1 ~ ^0 HV' 


k =0 


The generator truncated at order cc^ 4 " 1 is given by G[,v] = * ^ n=0 a n F™ +1 . 
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B. Computation of the ordered exponential 


In classical mechanics the ordered exponentials are usually computed using the Deprit 
“triangle” 19 . However, we have observed that for all of our quantum examples the following 
algorithm is faster. 


Up to O(o w+1 ) the transform of general perturbed Hamiltonian H = ct n H n can be 

written in the following form: 


N 


N 


H = exp_H / A? ds)H = nry v n \H = V n Ff). 

■'° \n =o / n =o 

Here we introduced ancillary double-indexed operators = J2k= n ak Hk-n- 

The Deprit relation (10) allows us to express repeatedly the superoperators V ra by means 

of ^less then n • 

N-l „ JV-1 ^ n 

H = vViv " 1 + V V*ff > = V V*Ff= ... = £ vTl 


T™) 


k =0 


fc =0 


fc =0 


where the operators F^\ n = N — 1,..., 0 are computed using the following relation: 


p( n )_ pD+i) 

r k ~ r k 


1 r p( n +!) 
' L 'G T ,,. k r n +1 > 


k — 0 ,..., n. 


n+1 n ~ k 

Finally, all the superoperators V disappear and we obtain the transformed Hamiltonian: 


— f(°) 


H = F, 


C. Examples and comparison with other methods 

The explicit expression is important from a general mathematical standpoint because 
it systematizes and simplifies the perturbation expansion. Moreover, it provides a suf¬ 
ficiently effective computational algorithm. Here we compare the computer times of 
block-diagonalization using the explicit expression (12) with those of the Van Vleck 4 and 
Magnus 27,28 approaches. 

The Van Vleck method block-diagonalizes the Hamiltonian using a chain of exponents: 

_ e ~ia n G n -i e ~iaG 0 e iaG 0 e ia n Gn-i _ e ~ to ‘ n ^G n _ 1 e ~ ia ^G 0 

In classical mechanics it is known as the Dragt-Finn transformation 29 . 
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The Magnus method uses the superoperator exponent with an Q-dependent generator 5 

fj = g—* Hi ak Gk~i jj g* H" ak Gk~i __ g-*E 


In classical mechanics this is known as the Hori algorithm 30 . 

We have evaluated the efficiency of these methods for the following model systems: 


• The one-dimensional oscillator with the quartic anharmonicity: 


H = 



The first orders of the block-diagonalized Hamiltonian are: 


H =h(N + —) + ah 2 


31 2 3 

-(N+-) H- 

8 V 2 32 


— a 2 h 3 


17, - l 3 67 , - 1, 

-(N+-) H- (N+- 

64 1 2 256 v 2' 


0(a 3 


Here N = a)a. This is the typical structure of the perturbation series for nondegenerate 
systems. It is worth noting that the effective diagonal Hamiltonians of such systems 
are identical for all the CPT methods. 


The Henon-Heiles system. This is a two-dimensional system with the Hamiltonian 

H = ^(Pi + ql+Pl + ql) + - ^ 2 )- 


The block-diagonalized effective Hamiltonian is: 


H = ( h--a 2 h 2 -—. 

9 108 


a 4 h 3 1 + I h — -a 2 h 2 — —a^h 


2fc2 


61 


4fc3 


54 


cqcq —I- a,a 2 


- ( ]h a2h 2 + S" 4;i3 ) + 

1 2»-2 9 4.3 \ 

-a n - a n alaia,a, 

3 4 11122 


al 2 d 2 2 | - 


7 2 1-2 7 4 3 

—ah - ah 6 

12 48 


~f2~2 1 ~t2~2\ 

6^2 CL-[ CL\ CL 2 ) 


+ a‘tf<7M 


,t3; 


65 


,t2 -1 A2; 


175 


a a?-a a,a,-afaia,a, + 

432 1 1 144 1 1 2 16 1 2 1 2 144 


a\ 2 d\dl 


235 


,t3~3 


161 ^f A t2 A 3 . 47 f A 2 1 175 ftt3ft2ft 1 1 


-a 9 “a5- ala, a? H- ala, a-,a, H-a, afa, + O(cr). 

432 2 2 144 1 2 1 16 1 2 1 2 144 2 1 2 J y ’ 

This expression could be further diagonalized using finite canonical transformations 31 . 


All three CPT methods build near-identity unitary block-diagonalizing transforma¬ 
tions. Such transformations are equivalent to ordered exponentials with different con¬ 
ditions 'PfjF. As a result, all the three effective block-diagonal Hamiltonians of the 
Henon-Heiles system differ from each other starting from the 8 th perturbation order, 
and are connected by block-diagonal unitary transformations. 
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The comparison of efficiency of these methods is not unambiguous. High order computa¬ 
tions process large multi-gigabyte expressions containing millions of terms. The performance 
of operations with such expressions depends on the computer algebra system used, its inter¬ 
nal optimizations, the server CPU, RAM, disk subsystem, etc. Even the relative efficiencies 
of methods may vary. Therefore, the following our results are only illustrative. 

Figure 1 compares the computational times of block-diagonalization for the Henon-Heiles 
system on an Oracle™ Exadata X2-2 server with Intel Xeon X5675 (3.06 GHz) processor 
using Form 4.1 computer algebra system 12 . 



The computation times for the quartic anharmonic oscillator follow the same pattern. In 
these particular cases the explicit algorithm is faster for high-orders computations than the 
Magnus expansion, but less efficient then the Van Vleck method. 

VI. SUMMARY 

We have presented here an application of Kato resolvent expansion to quantum superoper- 
atorial canonical perturbation theory. The canonical identity for the superresolvent allowed 
us to demonstrate unitary intertwining of perturbed and unperturbed averaging superpro¬ 
jectors. This leads to the explicit expression for the generator of the block-diagonalizing 
transformation in any perturbation order and systematic description of ambiguities in the 
generator and block-diagonalized Hamiltonian. 
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We have also compared the computational efficiency of the explicit expression for the 
quartic anharmonic oscillator and the Henon-Heiles system with the efficiencies of the Van 
Vleck and Magnus methods up to the 32 nd perturbation order. 
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